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Several  edge-on  impact  tests  on  transparent  glasses  and 
polycrystalline  ceramics  have  shown  that  failure  fronts  have  an 
extremely  rough  morphology,  including  appearance  of  spikes  and 
cracks.  A  simple  thermodynamic  theory  has  been  used  to  predict  the 
appearance  of  the  roughness  as  a  manifestation  of  morphological 
instability  of  failure  fronts.  For  the  case  of  isotropic  phases,  the 
instability  criterion  can  be  obtained  in  explicit  fonn. 


INTODUCTION 

When  a  high-speed  projectile  impacts  a  brittle  material,  such  as  glass  or 
polycrystalline  ceramic,  severe  damage  and  fragmentation  is  normally  observed  before 
projectile  penetrates.  AlON  is  a  material  being  considered  for  a  variety  of  transparent 
armor,  sensor  window,  and  radome  applications.  It  is  a  polycrystalline  ceramic  that 
fulfills  the  requirements  of  transparency  and  requisite  mechanical  properties  for 
transparent  armor  against  armor  piercing  ammunition.  AlON  has  a  cubic,  spinel  crystal 
structure  (Fd3m)  that  can  be  processed  to  transparency  in  a  polycrystalline 
microstructure.  It  differs  from  glasses  which  do  not  have  any  periodic  crystalline  order, 
but  it  is  akin  to  polycrystalline  opaque  ceramics,  such  as  aluminum  oxide. 

AlON  was  recently  investigated  by  Strassburger,  Patel,  McCauley,  and  Templeton 
(2006)  using  EOI  test  in  two  different  optical  configurations.  In  the  first,  a  regular 
transmitted-light  shadowgraph  set-up  was  used  to  observe  wave  and  damage 
propagation.  In  the  second,  a  modified  configuration  was  used,  where  the  specimens 
were  placed  between  crossed  polarizers,  and  the  photoelastic  effect  was  utilized  to 
visualize  the  stress  waves.  Pairs  of  impact  tests  at  approximately  equivalent  velocities 
were  carried  out  in  transmitted  plane  light  (shadowgraphs)  and  crossed  polarized  light. 
AlON  and  fused  silica  specimens  were  impacted  using  solid  cylinder  steel  projectiles 
with  velocities  ranging  from  270  to  925  m/s. 
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b)  cross -polarized  light 


Figure  1.  Selection  of  two  shadowgraphs  and  corresponding  crossed  polarizers 
Pictures  from  impact  on  AlON  at  380  m/s. 

The  nucleation  of  crack  centers  was  observed  ahead  of  the  apparent  fracture  front, 
growing  from  the  impacted  edge  of  the  specimens.  A  comparison  of  the  shadowgraphs 
the  photographs  (recorded  in  a  reflected  light  configuration  with  a  coated  AlON 
specimen  at  the  same  impact  conditions)  indicated  fracture  nucleation  in  the  interior  of 
the  ceramic. 

Figure  1  shows  a  selection  of  two  regular  shadowgraphs  along  with  the 
corresponding  crossed  polarizer  photographs  for  tests  with  AlON  at  380  m/s.  The  high¬ 
speed  photographs  show  rapidly  growing  darkened  to  opaque  regions,  which  reflect 
changes  in  the  optical  transmission  due  to  pressure-induced  refractive  index  changes 
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and  damage  and  fractured  zones  within  the  specimen.  In  addition,  the  nucleation  of 
crack  centers  ahead  of  the  crack  front  is  clearly  visible  8.7  ps  after  impact.  In  contrast  to 
the  shadowgraphs,  where  a  wave  front  is  not  discernible,  the  crossed  polarizer’s 
configuration  reveals  an  approximately  semicircular  stress  wave  front,  which  is  a  little 
further  advanced  than  the  damage  front  visible  in  the  shadowgraphs  at  the  same  time.  In 
the  shadowgraph  images  (Figure  la)  the  failure  (damage)  region  clearly  exhibits  very 
sharp  failure  and  damage  spikes  extending  out  from  the  front  of  the  main  damage 
region. 

2.  FORMAL  STATEMENT  OF  THE 
PROBLEMEquation  Chapter  (Next)  Section  2 

In  the  paper,  we  suggest  a  simple  thermodynamic  theory  which  explains  the 
appearance  of  the  spikes  in  the  aforementioned  experiments  with  AlON.  Two  different 
theoretical  and  computational  approaches  seem  to  be  most  useful  in  analyzing  different 
aspects  of  wave  propagation  in  brittle  solids.  They  describe  failure  fronts  on  different 
spatial  length  scales.  According  to  the  first,  less  detailed,  “macroscopic”  approach,  the 
failure  front  is  treated  as  an  infinitely  thin  mathematical  surface  separating  intact  and 
comminuted  spatial  domains.  According  to  the  second,  “microscopic”  approach,  the 
failure  front  has  very  thin  (though  still  finite)  thickness.  Within  this  thin  spatial  domain, 
all  physical  parameters  of  the  substance  change  continuously  but  extremely  rapidly. 
Quite  often,  additional  physical  processes  should  be  taken  into  account  inside  these  thin 
domains,  whereas  outside  these  thin  domains,  the  additional  physical  processes  play  a 
relatively  minor  role.  In  this  respect,  the  situation  with  the  two-scale  description  of 
failure  waves  is  the  same  as  with  the  classical  hydrodynamic  and  combustion/ 
detonation  shock  waves  (Landau  and  Lifshitz.,  1987;  Courant  and  Friedrichs,  1948). 

Although  the  spatial  domain  of  intensive  fracture  in  the  Edge-on  Impact  AlON 
experiments  are,  in  fact,  quite  thick,  we  begin  our  analysis  with  the  simpler  macroscopic 
approach.  Within  each  of  the  spatial  domains,  the  energy  densities  eml  and  edam  per 
unit  mass  of  the  intact  and  damaged  states  are  given  by  the  following  quadratic  forms  of 
displacement  gradients,  respectively: 


Pein{umn)  =  {\l2)c 


ijkl 

■  ,  U  U,  , , 
int  i,j  k,l  ’ 


PedaJUm,n)  =  (]/2)c 


famUiJUk,l  +  Clb 


(2.1) 


where  cf*7  and  cf’m  are  the  elasticity  tensors  of  the  two  states,  and  p  is  the  original 
mass  density  (a  subscript  after  comma  is  used  to  denote  differentiation  with  respect  to 
the  spatial  coordinates  x' ).  Limiting  ourselves  with  the  approximation  of  linear 
elasticity,  we  ignore  all  effects  of  mass  density  change.  The  positive  constant  qh  takes 
into  account  the  energy  required  to  produce  various  defects  (interfaces,  vacancies,  shear 
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bands,  holes,  etc.)  distributed  within  unit  mass  of  the  bulk  of  damaged  substance.  This 
term  is  analogous  to  the  constant  used  in  the  theory  of  slow  combustion,  which  takes 
into  account  the  energy  release/consumption  due  to  chemical  reactions  (Landau  and 
Lifshitz.,  1987;  Courant  and  Friedrichs,  1948)..  See  the  papers  (Grinfeld,  and  Wright, 
2002,  2004)  for  a  more  detailed  discussion  of  the  model. 

In  addition  to  appropriate  initial  and  external  boundary  conditions,  the  master 
system  includes: 

a)  the  bulk  momentum  equation  within  each  of  the  bulk  domains, 


d2u  _  dpJl 


(2.2) 


b)  the  displacement  continuity  equation  across  the  failure  front, 

[«']-«,  =  0  , 


(2.3) 


c)  and  the  momentum  continuity  across  the  failure  front, 


(2.4) 


where  pJl  =  de(umn)/ 8uUj  is  the  stress  tensor,  and  ni  is  the  unit  normal  to  the  failure 
front. 

The  last  equation  across  the  failure  front  describes  the  kinetics  of  failure  as  follows: 


c 


(2.5) 


where  c  is  the  velocity  of  the  failure  front,  and  the  tensor  p  jk  is  defined  as 

pjk  =  eS[  -  p  pji  (Sik  +  «,  )  .  (2.6) 

It  plays  the  same  role  as  the  scalar  Gibbs  chemical  potential  p  of  a  liquid  substance; 
K  is  a  positive  (kinetic)  constant  or  function  with  dimension  [velocity]  '.  We  assume 
explicitly  that  the  displacements  and  traction  remain  continuous  across  the  interface. 
This  is  not  the  only  reasonable  option.  Another  reasonable  option,  especially  when 
dealing  with  highly  fragmented  states,  would  be  the  model  of  a  friction-free  interface 
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with  discontinuous  displacements.  In  this  case  the  last  (kinetic)  constitutive  equation 
(2.5)  should  be  modified  as  well  (a  similar  system  was  analyzed  in  the  context  of  phase 
transformations  in  the  monograph  (Grinfeld,  1991)). 

We  limit  our  study  with  two-dimensional  propagation  and  consider  an  initially 
resting  uniform  half-plane,  x  >  0 ,  experiencing  an  impact  at  x  =  0  by  the  oblique  force 
P  (see  Figure  2). 


z 


Damaged  1 
substance  ' 

Intact 

L  substance 

V, 

/ 

|  Failure  front 

Figure  2.  Oblique  impact  of  a  brittle  substance. 


The  system  (2.2)-(2.6)  allows  piece-wise  linear  solutions  of  the  following  form: 

u[(x,z,t)  =  d\t  +  a\z  at  x>  ct ,  u'_(x,z,t)  =  d't  +  a_z  at  x  <  ct ,  (2.7) 


where  d‘±,  a\  ,  and  c  are  certain  constants.  The  sign  +  (-)  marks  the  quantities  related 
to  the  intact  (damaged)  state.  The  constant  c  gives  the  speed  of  the  failure  front  across 
which  the  displacement  gradients  a'±  suffer  the  finite  discontinuities.  The  solution 
mentioned  above  allows  the  consideration  of  the  quasi-static  problem  oblique  impact 
(loading)  when  the  applied  force  has  the  components  p±  =  n  and  p+  =  r  .  The 
resulting  formula  is  particularly  simple  and  transparent  in  the  quasi-static 
approximation,  i.e.,  when  inertia  is  neglected  in  the  equations  (2.2)  and  (2.4);  then, 
formula  of  the  velocity  of  the  failure  front: 


n1 

r  1 

1  ] 

T2 

_l _ 

n  _n 

2 

(  A_  +2  ju 

+  2p+  J 

i 

2 

(2.8) 


here  X±  ,  /u±  are  the  Lame  elastic  modules  of  the  intact  and  comminuted  materials  that 
are  assumed  isotropic.  The  Poisson  ratios  are  defined  as  v±  =  X±  /  2(/t±  +  //±  ). 
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3.  MORPHOLOGICAL  INSTABILITY  OF  FAILURE  FRONT 
Equation  Section  (Next) 

In  order  to  explore  morphological  stability  of  the  piece-wise  linear  solution  we 
present  the  elastic  displacement  u‘±(x,z,t )  and  the  speed  of  the  interface  c(z,t)  in  the 
following  fonn: 

u‘±(x,z,t)  =  u+  (x,t)  +  ul(x,z,t),  c(x,z,t)  =  c°  +c(z,t),  (3.1) 

where  u'±(x,z,t )  and  c(z,t)  are  small  disturbances.  The  above  presentations  should  be 
substituted  in  the  bulk  equations  and  in  the  boundary  conditions  (1),  (2),  which  should 
then  be  linearized  with  respect  to  the  small  disturbances.  We  then  look  for  the  solutions 
of  the  linearized  system  in  the  following  form  (k  is  the  in-plane  wave-number  and  //  is 
the  rate  of  growth): 

(3.2) 


u 


I (x, z, t)  =  W‘  (jc  -  c°t)eik+* ,  c(z, t)  =  Seikz+tn 


where  the  functions  W+  ( X )  should  decay  exponentially  at  X  — »  ±  oo . 

The  dispersion  equation  for  r)  is  somewhat  lengthy  and  cumbersome  without  further 
special  assumptions.  For  the  case  of  isotropic  states  and  two-dimensional  geometry  it 
reads  as  follows: 


K  4//  (//2  -//2)(l-u  ) 

k  n  ~  (3  -  4v_ )  K  +  Tf  A  (1  -  v  )  +  K 
_ (1  -  v_ )  +  JU- _ 

(3-4v_)/ul  +4fi_fi+(l-v_)  +  ti*  ' 


-8 


r2  n_{n_  -//+)2  (1  —  )2 

M+  (l-2u)2 


(3.3) 


The  last  formula  becomes  relatively  compact  and  instructive  in  the  limit  of  an 
incompressible  damaged  state  (v  =  1  /  2  )  and  when  in  the  intact  state  v+  =  1  /  3  .  Then, 
it  reads  as  follows: 

M+  (/L2  -M-)  1  r 2  /L  ju+  ~ )~  4//  +  3 /u+ 


P 

8  Kk 


1 1  = 


5ju_  +Sju_ju+  +3ju+  n 


P 


5jU2  +  8  jU_jU+  +  3  jui 


(3-4) 


The  wave  front  is  morphologically  unstable  if  the  linearized  master  equation  permits 
solutions  of  the  above  form  and  with  the  rate  rj  corresponding  to  exponential  growth  in 
time.  The  last  fonnula  (3.4)  shows  that  shear  stresses,  r  ,  play  a  stabilizing  role  (similar 
to  the  case  of  morphological  instabilities  of  solid-solid  phase  interfaces,  discussed  in 
Grinfeld,  1991).  In  addition,  at  r  =  0,  the  failure  front  is  morphologically  unstable  in 
the  most  interesting  case  when  the  shear  modulus  /u  of  the  damaged  state  is  less  than 
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the  shear  modulus  ft+  of  the  intact  state.  This  instability  has  a  simple  physical  meaning. 

It  means  that  penetration  of  fingers  of  damaged  material  into  intact  material  is  the 
fastest  way  of  releasing  accumulated  elastic  energy  from  the  system. 


4.  INTUITIVE  INTERPRETATION  OF  THE  MORPHOLOGICAL 
INSTABILITY  OF  FAILURE  FRONTS 

A  plausible  interpretation  of  the  stress  driven  destabilization  of  failure  fronts  can  be 
based  on  the  principle  of  the  fastest  energy  release  (or  of  the  entropy  production)  of 
irreversible  thennodynamics.  In  the  context  of  fracture,  the  principle  postulates  that 
fracture  develops  in  such  a  way  that  the  maximum  elastic  energy  should  be  released  (per 
fixed  amount  of  damaged  bonds).  Simply  stated,  the  propagation  of  fingers  has  a  simple 
physical  meaning.  It  means  that  penetration  of  fingers  of  damaged  material  into  intact 
material  is  the  fastest  way  of  releasing  accumulated  elastic  energy  from  the  system. 


a)  b) 

Figure  3.  Towards  intuitive  interpretation  of  the  hackle  zone  formation. 

Figure  3  shows  two  different  morphologies  of  the  advancing  front:  Figure  3  a 
illustrates  an  advancing  failure  front  that  does  not  change  its  morphology  and  remains 
flat,  and  the  Figure  3b  shows  an  advancing  front  in  the  form  of  a  horizontal  finger.  It  is 
assumed  that  the  two  newly  created  comminuted  areas,  shown  in  red  in  the  figures  3  a 
and  3b,  are  equal.  It  should  be  intuitively  clear  that  the  finger-like  morphology  of  the 
failure  front  allows  much  more  elastic  energy  to  be  released  than  the  morphologically 
stable  front  shown  on  3a).  Indeed,  the  propagating  flat  front  on  3a)  releases  the  elastic 
energy  only  within  the  small  rectangle  (between  the  original  and  final  front  positions), 
whereas  the  finger-like  front  allows  to  release  the  accumulated  elastic  energy  from 
everywhere  to  the  right  of  the  original  front  (including  in  the  material  that  remains  in 
the  intact  state).  So  the  "fingering"  of  the  front  allows  the  accumulated  energy  to  be 
released  in  a  more  efficient  way  than  with  the  morphologically  stable  propagation  of  the 
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flat  failure  front.  This  is  the  potentially  possible  mechanism  of  morphological 
destabilization  of  the  failure  front  discussed  in  this  paper. 

5.  CONTINUUM  MODEL  OF  DAMAGE 
Equation  Section  5 

The  model  considered  above  is  a  certain  two-state  discrete  fracture  model  (Grinfeld 
and  Wright,  2002,  2004).  A  continuum  fracture  model  should  include  a  damage 
parameter  k  that  can  assume  a  continuum  range  of  values.  Corresponding  elastic 
energy  e  function  is  becoming  a  function  of  the  elastic  gradients  nl|y  and  of  the  damage 

parameter  k:  e  =  e(uij,/c).  For  such  a  model  the  bulk  equations  of  equilibrium 

includes  a)  the  condition  of  mechanical  equilibrium,  and  b)  the  equation  of  “chemical” 
equilibrium,  i.e.,  of  equilibrium  with  respect  of  damage  proliferation: 


a) 

b) 


8  de{umn,K) 
dxJ  duij 

,„,*•)  _Q 

8k 


(5.1) 

(5.2) 


It  is  more  or  less  typical  for  the  system  with  damage  that  full  equilibrium  with 
respect  to  damage  proliferation  cannot  be  achieved.  In  such  situations,  we  have  to 
switch  from  the  equation  of  chemical  equilibrium  (5.2)  to  its  quasi-static  analogue: 


8k  8e{um  k) 

—  =  —T) - 

8t  8k 


(5.3) 


where  rj  is  a  certain  positive  kinetic  constant. 


The  important  issue  here  is  the  choice  of  an  appropriate  energy  function 
e  =  e(uiU,K) .  For  the  purposes  of  fracture  mechanics  such  a  function  can  be  found  by 

combining  the  approaches  of  the  theories  of  vacancies  in  ceramics  sintering  (Lifshitz, 
1963)  and  of  the  damage  mechanics  (Kachanov,  1986). 

The  quasi-static  system  presented  above  was  used  to  numerically  explore  the 
behavior  of  a  circular  plate  with  a  thin  elliptic  cavity  under  action  of  external  pressure 
for  the  energy  density 
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e(uiU,/c)  =  ^(/c-K°)2  +  ju</){k) 


- U  UJ  .  +U,,  ;.U 

\-2v  '•  J  °’7)  ' 


(5.4) 


with  the  damage  function  (fftc)  (Grinfeld  and  Wright,  2002,  2004) 


=  1  -  (l  -  cmin  )-^r ,  0  <k<k,  0  <  cmin  <  1 . 

K 


(5.5) 


Thus,  the  suggested  model  depends  on  6  constants:  //,  v,  q,  k  ,  k*  ,  cmin.  The  first 
two  of  them  are  just  the  shear  modulus  and  Poisson’s  ratio  of  the  intact  substance.  The 
physical  meaning  of  the  remaining  four  constants  will  be  explained  elsewhere.  For  a 
certain  choice  of  the  constants,  the  evolution  of  the  relative  damage  parameter  k  l  k  *  is 
shown  in  the  Figure  4.  The  depth  of  gray  color  reflects  the  magnitude  of  the  relative 
damage  parameter  —  white  corresponds  to  a  fully  damaged  material  and  black  to  a  fully 
undamaged.  Not  surprisingly,  the  largest  damage  occurs  at  the  tips  of  the  cavity. 


Figure  4.  Development  of  damage  in  a  circular  elastic  plate  with  an  elliptic  cavity. 

Figure  5  illustrates  developing  damages  for  the  case  of  a  brittle  ring  under  the  action 
of  different  outer  and  inner  pressure.  It  is  very  instructive  from  at  least  two  points  of 
view. 

First,  this  geometry  is  representative  of  the  important  problem  of  penetration  of  a 
rod- like  projectile  with  a  circular  cross  section  through  a  thin  brittle  isotropic  ring. 

Second,  this  example  is  particularly  instructive  from  the  standpoint  of 
morphological  instability  of  failure  fronts.  The  formulated  problem  is  obviously  radially 
symmetric  since  the  intact  and  damaged  states  are  assumed  isotropic,  the  boundary 
circles  are  assumed  concentric,  and  the  pressure  loading  does  not  violate  the  radial 
symmetry.  Thus,  the  problem  under  consideration  possesses  a  solution  describing 
evolution  with  radial  symmetry.  This  solution  would  be  easy  to  model  numerically  if  it 
were  stable  with  respect  to  small  perturbations  of  radial  symmetry,  however,  such  non- 
radial  small  perturbations  are  always  present  in  any  numerical  implementation. 
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Figure  5.  Isotropic  brittle  ring  under  unequal  outer  and  inner  pressures. 

The  results  of  (Grinfeld  et  al,  2006),  briefly  presented  in  the  section  3,  prompt  us  to 
make  a  conjecture  that  the  radial  solution  of  the  system  (5.1),  (5.2)  for  the  ring  under 
unequal  pressures  is  unstable  and  cannot  be  obtained  numerically.  Figure  6  completely 
confirms  this  conjecture.  At  this  early  stage  of  our  research,  it  remains  totally  unclear 
what  parameters  control  the  number  of  microcracks  in  the  model.  It  will  be  the  subject 
of  our  future  theoretical  analysis  and  numerical  experiments. 

6.  CONCLUSION 

A  transparent  polycrystalline  AlON  was  tested  in  an  Edge-On  Impact  configuration. 
Damage  propagation  was  observed  in  a  regular  shadowgraph  set-up  and  in  a  dynamic 
photoelastic  model  using  crossed  polarizers.  Photographs  recorded  in  a  reflected  light 
configuration  with  a  coated  AlON  specimen  at  the  same  impact  conditions  confirmed 
that  damage/fracture,  and  isolated  crack  nucleation  occurred  in  the  interior  of  the 
ceramic;  the  failure  (damage)  front  proceeded  through  the  AlON  plate  with  a  spikes-like 
morphology. 

We  have  presented  a  simple  model  for  a  failure  front  based  on  an  analogy  with  slow 
combustion,  and  we  have  presented  a  formula  for  the  velocity  of  failure  wave  generated 
by  oblique  impact  on  a  brittle  material.  We  then  demonstrated  that  under  rather  general 
assumptions,  the  flat  failure  front  is  morphologically  unstable.  The  last  conclusion 
shows  an  additional  analogy  between  failure  fronts  and  slow  combustion  fronts 
(however,  the  mechanisms  of  destabilization  are  totally  different). 

Theoretical  results,  based  on  the  two-state  model  of  damage,  have  been  confirmed 
by  numerical  experiments  based  of  the  IDZ-models  suggested  earlier  in  (Grinfeld  and 
Wright,  2002,  2004) 
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The  morphological  instability  of  failure  fronts  sheds  light  on  a)  the  appearance  of 
spikes  in  the  AlON  edge-on  impact  experiments,  b)  on  failure  fronts  and  crack 
bifurcations  as  observed  in  experiments  on  static  indentation  of  brittle  materials,  and  c) 
on  appearance  of  radial  cracks  in  dynamic  experiments  with  penetration  of  projectiles 
through  brittle  materials. 

This  suggests  that  the  investigation  into  AlON  offers  insight  into  damage  evolution 
not  only  in  AlON,  but  also  in  other  opaque  ceramics  where  only  surface  damage  can  be 
observed  with  optical  methods.  Further  experimental  and  theoretical  studies  are 
mandatory  before  making  ultimate  conclusions. 
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